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Baryons in Partially-Quenched Chiral Perturbation Theory 

M. J. Savage a 

department of Physics, University of Washington, Seattle, WA 98195-1560, USA 

I discuss the inclusion of baryons into partially-quenched chiral perturbation theory and describe one-loop 
calculations that have been performed. 



1. Introduction 

The motivation for including baryons into 
partially-quenched QCD (PQQCD), or more 
specifically partially-quenched chiral perturba- 
tion theory (PQxPT), is clear. One wishes to 
take partially-quenched (PQ) lattice simulations 
of observables of interest that in the foreseeable 
future will be performed with quark masses that 
are heavier than those of nature, and make rig- 
orous predictions for these observables in QCD. 
The way that such predictions will be made is to 
use PQ lattice simulations to "match" onto a PQ 
effective field theory (PQEFT), and then use the 
explicit quark mass dependence of the PQEFT 
to make predictions at the physical values of the 
light quark masses, i.e. the coefficients of the lo- 
cal operators in the PQEFT will be determined 
from the PQ lattice simulations. 

1.1. Quenching and Partial- Quenching 

Ideally, all lattice simulation would be un- 
quenched and the quarks would have their phys- 
ical masses. At some point in the future this 
will happen, but at this point in time and for 
the foreseeable future such simulations will not 
be possible, simply due to the lack of comput- 
ing power. It is the quark loops that are dis- 
connected from external sources (the quark de- 
terminant), participating through gluons alone, 
that are the most time consuming to simulate as 
they scale as a relatively high power of the inverse 
quark mass. In contrast, simulation of quark 
loops that are coupled to the external sources is 
much less time consuming. In quenched QCD 
(QQCD) |l|-Q], the quark determinant is set equal 



to unity, and thus the time-expensive discon- 
nected quark loops are discarded, allowing for rel- 
atively rapid simulation at relatively light quark 
masses. From a formal standpoint, the quenched 
simulations correspond to introducing an equal 
mass "ghost" fermion for each light quark, so that 
disconnected loop contributions from the light 
quarks and the ghosts exactly cancel. Thus the 
SU(2) L (g> SU(2) R chiral symmetry of QCD is ex- 
tended to SU(2\2) L ® SU(2\2) R x 17(1). There 
is no formal limit in which quenched simulations 
can be used to make rigorous predictions about 
QCD Q The discovery of PQQCD [|§ was a 
major step forward in this regard. In addition 
to the ghosts that are added in QQCD, addi- 
tional quarks are added to partially unquench 
the simulation, increasing the flavor symmetry 
of the theory to SU(A\2) L <g> SU(A\2) R (g> U(l). 
These additional "sea" -quarks participate in the 
disconnected loop diagrams, but if they are heav- 
ier than the "valence" -quarks they require less 
time to simulate. In the limit that the mass 
of the sea-quarks become equal to those of the 
valence-quarks (and ghosts) the matrix elements 
of PQQCD between physical states become equal 
to those of QCD. Thus one can make QCD pre- 
dictions from PQQCD. 

In order to describe the low-momentum dy- 
namics and the light-quark mass dependence of 
the low-lying hadrons in QCD one, of course, uses 
an EFT, such as xPT in the meson sector and 
heavy baryon %PT (HBxPT) in the single baryon 
sector. In extending the theory to PQQCD one 



however, it has recently been noted that QQCD in 
the large-TVc limit is identical to QCD in the large-TVc 
limit g|. 
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has to extend xPT to PQxPT. This allow one to 
determine the contributions that are non-analytic 
in external momentum and m q , and also to pa- 
rameterize short-distance physics in the local op- 
erators that enter order-by-order in the expan- 
sion. As it is presently quite difficult to isolate 
the non-analytic contributions from the analytic 
ones in lattice simulations due to the relatively 
large values of m q that can be simulated, knowl- 
edge of the low-energy EFT is vital. The lattice 
simulations are used to determine the values of 
counterterms that appear at any given order in 
the EFT. Thus, PQ simulations will be used to fix 
counterterms in PQxPT, which can then be used, 
to make QCD predictions as the m g -dependence 
is explicit. Such determinations are beginning to 
be performed in the meson sector 

2. Baryons in PQxPT 

Baryon properties in quenched %PT (QxPT) 
have been studied somewhat during the past 
six years, starting with the seminal paper by 
Labrenz and Sharpe [Q. The representations of 
the baryons are extended from those of ordinary 
lie-groups to those of graded lie-groups, SU(3) 
to S[/(3|3) for three-flavor flavor QCD and from 
SU(2) to SU(2\2) for two-flavor QCD, so that the 
tools developed for HBxPT can be used directly 
in PQHBxPT. In the quenched theory, the octet- 
baryon masses [|J , magnetic moments Q and the 
matrix elements of isovector twist- 2 operators j^] 
have been studied at the one-loop level in the chi- 
ral expansion. 

Including baryons in PQxPT is similar to in- 
cluding them into Q%PT pO| , pl| . One constructs 
the irreducible representations of the graded lie- 
group, which for two-flavors is SU(A\2), and as- 
signs baryon fields, just as one does in QCD 
and QQCD. The 70 dimensional representation 
of baryons that contains the nucleons is described 
by a tensor of the form Bijk , where i, j, k run from 
1 to 6 (we have dropped the Dirac index) . To de- 
termine the particle assignment it is particularly 
simple to use the interpolating fields discussed by 
Labrenz and Sharpe, 



for the representation containing the nucleons. 
The objects Qi are quark field operators that in- 
clude the valence, sea and ghost quarks. C is 
the charge conjugation operator, a, b, c are color 
indices and a, f3, 7 are Dirac indices. The interpo- 
lating field for the 44 dimensional representation 
that contains the 3/2-baryon resonances, Tijk (we 
have dropped lorentz and Dirac indices), is 



QT a Qf b QV + cyclic e abc (cr)^ 



(2) 



QV 



7.bc (C75) a/3 (1) 



One finds that under the interchange of flavor in- 
dices m, 

Bi jk = R 1+ "^ B ikj (3) 
Ban + R 1+ '™ Bjik + 

Hjk = (-) 1+ViVj T jik = {-Y+wH T . k . . 

with r/k — +1 for k = 1,2,3,4 and 77^ = for 
k = 5,6. 

The only baryons that appear as intermediate 
states in one-loop diagrams involve either three 
valence quarks, or two valence quarks and one 
ghost quark, or two valence quarks and one sea 
quark. Thus we only need to determine the lo- 
cation of these baryons in the tensors B and 
T. It turns out to be convenient to classify the 
baryons according to how they transform under 
SV(2)vaicncc ® SU(2) sca ® S'J7(2) ghost , and it is 
straightforward to show that for the 70 one needs 
to include fields N a , a t, a Sbc, at an d a Sbc that 
transform as (2,1,1), (1,2,1), (3,2,1), (1,1,2) 
and (3,1,2) respectively. For the 44, one intro- 
duces fields A a fo c , a Xbc and a Xb c that transform 
as (4,1,1), (3,2,1) and (3,1,2) respectively. For 
the three-flavor case the octet of baryons are part 
of a 240 dimensional irrep of SU(6\3) while the 
decuplet of baryon resonances are part of a 138 
dimensional irrep. 

It is straightforward to construct the lagrange 
density describing the baryons and their in- 
teractions with the pseudo-Goldstone bosons 
associated with the spontaneous breaking of 
SU(A\2) L <g> SU(A\2) R down to SU(A\2) V , Thus 
the loop expansion for the low-energy properties 
of the nucleons in PQxPT is the same as that in 
xPT except that the number of particles that can 
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participate in loops is much larger, and one needs 
to keep track of minus signs associated with the 
loops involving fermionic mesons. 

Calculations of several observables at the one- 
loop level, both in two- and three-flavor PQxPT, 
have recently been completed, in addition to a 
tree-level analysis of the nucleon-nucleon (NN) 
potential. A summary of what presently exists 
can be seen in Table [jj In the two-flavor calcula- 
tions, both the valence- and sea-quarks are non- 
degenerate, while in the three-flavor calculations 
the two light-quarks are taken to be degenerate 
(both valence and sea) while the strange quark is 
non-degenerate in both, i.e. 2 + 1. 

Instead of detailing each of these calculations, 
I wish to focus on a few aspects and observables 
that are somewhat out of the ordinary. For fur- 
ther details of the mainstream observables, such 
as the masses, magnetic moments and so forth I 
refer the reader to the appropriate papers. 

2.1. The Nucleon-Nucleon Potential 

In QCD it is well established that one pion ex- 
change (OPE) gives the dominant contribution 
to the long-distance component of the NN (NN) 
potential. OPE yields a NN potential that is 
Yukawa potential at large-distances, falling like 
~ e~ m * r /r and also provides a ~ 1/r 3 short- 
distance tensor contribution that mixes S-wave 
and I?- waves. 

Efforts have been made to explore the interac- 
tions between nucleons on the lattice by deter- 
mining their scattering lengths. However, this is 
exceptionally challenging due to the fact that the 
scattering lengths are unnaturally large, ~ 6 fm 
in the 3 Si — 3 Di coupled channels and ~ 24 fm 
in the ^So channel. What's more the simulations 
have been performed in QQCD |Q and currently 
no unquenched or PQ calculations exist (for a re- 
cent discussion see Ref. Q). 

In both QQCD and PQQCD the long-distance 
component of the NN potential is very different to 
that in QCD. The hairpin interaction that gives 
rise to a double pole in the iso-singlet propaga- 
tor, that in QQCD is proportional to Mq, and in 
PQQCD is proportional to m 2 v — m 2 s and van- 
ishes in the QCD limit, provides the dominant 
contribution to the long-distance component of 



the NN potential [Q. In the isospin limit, the 
SU{2) singlet (77, distinct from the 5Z7(4|2) sin- 
glet) propagator has the form 



i(m 2 sa - m 2 vv ) 



(4) 



{q 2 — m 2 v + if) 2 

where v = u = d denotes valence-quarks and 
s = j = I denotes sea-quarks. The mass m vv de- 
notes the mass of a meson containing two valence- 
quarks, and m sv denotes the mass of a meson 
containing one valence-quark and one sea-quark. 
Fourier transforming the momentum-space NN 
potential gives 



y(PQ)^ = °\ iXfg ' V e -"V 



STTf 2 



(5) 



2 T l ' T 2 

9 A 



9l 



(m 2 . - ml) 



277J„ 



which does not exhibit Yukawa type fall off 
at large distance but rather falls exponentially. 
Therefore, unfortunately, efforts to observe the 
long-distance behavior of the NN interaction in 
QQCD and PQQCD will need to subtract this 
large contribution. 

2.2. Charges 

When one considers processes involving elec- 
troweak gauge fields, there is an additional free- 
dom in charge assignments that is not present 
in QCD. One requires that QCD is recovered in 
the limit where the mass of the sea-quarks be- 
come identical to the mass of the valence quarks, 
and thus the electromagnetic charge matrix in 
QCD, QQ CD = diag(+§,-§), is extended to 
Q = diag(+§,-!,g,-,ffj,fc-,gj) in PQQCD. The 
charges qj, qi are arbitrary, and in the QCD 
limit, observables become independent of them. 
However, away from the QCD limit quantities 
do depend upon them, and it may be useful to 
make specific choices for particular observables. 
For instance, choosing qj = qi = means that 
disconnected loops coupling to a photon involve 
only valence-quarks, while choosing qj = +| and 
qi = — g means that disconnected loops coupling 
to a photon involve only sea-quarks. Thus, nu- 
merically, there may be some advantage in choos- 
ing the latter set of charges over the former, but 
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Table 1 



One-loop calculations in PQxPT that currently exist. 



Observable 


N f = 2 


N f = 3 


Masses 


V 


V 


Magnetic Moments @0 






Matrix Elements of Twist-2 Operators @@ 






Matrix Elements of the Axial Current JllJ 




X 


NNtt Parity- Violating Interaction 




X 


Nucleon Anapole Moment Jl(| 


V 


X 



this remains to be explored. In addition, there 
are sets of charges that minimize the contribution 
from one-loop diagrams in PQ%PT. As higher or- 
ders are parametrically suppressed, this will lead 
to the smallest variation in the matrix element 
with respect to changes in m q . 

The non-uniqueness in the extension of elec- 
troweak operators from QCD to PQQCD is not 
confined to the electric charges. It is also present 
for the twist-2 operators M], and for four-quark 
operators, such as those responsible for K — > 
7T7T [^5) and for NNtt parity- violation Jig . 

2.3. Proton Magnetic Moment 

As an example of a calculation in PQxPT, we 
discuss the proton magnetic moment, up to 

1 /2 

order 0{vhq ) in the chiral expansion. The lead- 
ing order contributions to the nucleon magnetic 
moment arise from the dimension-5 operators 

£(5) = 4^ F » u ^ Q (Ba^BQ) (6) 
+ (Ba^QB) + M 7 (Ba^B) strfi] , 

where the brackets (..) denote contractions of the 
flavor indices as discussed in Rcf. Q, and "str" 
denotes a supertrace. At tree-level, the coeffi- 
cients (J-ayfJ-p and /i 7 are related to the isoscalar, 
/j,o, and isovector, fii, magnetic moments of the 
nucleon, 



Mo = ~ {Ha + M/3 + 2 M 7 ) 
/ii = i (2^ a - fj,p) . 



(7) 



1/2 

Up to 0(rriQ ), the magnetic moment of the pro- 
ton can be written as 

M N 



where 



a P = Mo + Mi 

0p 



(8) 



(9) 



g\m vv 



[m s 



711 vv j 



\9a 



9 



9A9i 
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+qji (m 



2 2 1 1 2 

o9a + o9Agi + t9\ 
6 6 6 



0'p — ~^9an 



^•F vv Qjl 2 vv 3~ Sv) 



with 



11- 



The function Ti, 



T(mij , A, [i) is given by 



(10) 



where r\ = V A 2 — m 2 + ze, A is the A-nucleon 
mass splitting in the chiral limit, and \x is the 
dimensional regularization renormalization scale. 
qa is the usual axial coupling constant, while g\ 
is an additional axial coupling that contributes in 
PQxPT, but whose contribution must vanish in 
the QCD limit. gAN is the AiV7r coupling con- 
stant. 

There are two distinct contributions from the 
loop diagrams away from the QCD limit. One is 
a contribution that is independent of the choice 
of charges that vanishes like m sv — m VVl while 
the other also vanishes like m sv — m vv but de- 
pends upon the choice of charges of the sea- and 
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ghost-quarks. Therefore, one is able to find sets 
of charges for which the one-loop contribution to 
the proton magnetic moment from intermediate 
states in the 70 vanish, however, due to the non- 
trivial mass-dependences arising from intermedi- 
ate states in the 44, the one-loop contributions 
cannot be eliminated entirely. 

2.4. NNtt Parity Violation 

While flavor-changing parity- violating (PV) in- 
teractions are well understood theoretically and 
a great deal of precise data exists, knowledge of 
flavor-conserving parity-violation is rather primi- 
tive. Flavor-conserving parity-violation continues 
to be an area of intense investigation in the nu- 
clear physics community. Its study is presently 
serving both to uncover the structure of the 
nucleon in electron-scattering experiments such 
as SAMPLE 0, and to determine PV flavor- 
conserving couplings between pions and nucle- 
ons ||l8| , |l9| . The problems that are encountered 
in this sector are both experimental and theoret- 
ical. On the experimental side, the PV signals, 
unlike those in flavor-changing processes, appear 
as small deviations in either a strong or an elec- 
tromagnetic process, such as PV in ep — > ep, 
or in the circular polarization of the 7-ray emit- 



ted 
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The current situa- 



tion is somewhat confused by the fact that mea- 
surements of parity-violation in atoms and nu- 
clei do not give rise to a consistent set of cou- 
plings between hadrons pOfl . However, it is im- 
portant to keep in mind that many of the "experi- 
mental" determinations of these couplings require 
theoretical inputs with varying degrees of relia- 
bility. Recently, it has been reemphasized that 
measurements of PV observables in the single- 
nucleon sector would significantly ameliorate the 
situation by eliminating many-body uncertain- 
ties |^l|,^2| • Despite the inherent difficulty of such 
experiments, there are ongoing efforts to mea- 
sure PV processes in systems with only one or 
two nucleons, such as the angular-asymmetry in 
ftp — > e?7 p3| . Such measurements should pro- 
vide a reliable determination of the leading-order 
(LO), momentum-independent weak irNN cou- 
pling constant, h£jf N . 

On the theoretical side, despite heroic efforts 



to model [|2j,f25| hadronic matrix elements of 
the four-quark operators that appear in the low- 
energy effective theory of the standard model, 
there are no reliable calculations of the PV cou- 
plings between hadrons. A first principles calcu- 
lation of /i^jvjv m lattice QCD would therefore 
be extremely welcome. This would require a lat- 
tice QCD simulation of a correlator with three 
hadronic sources interacting via a four-quark op- 
erator. Unfortunately, chiral symmetry does not 
allow one to relate the wNN correlator to a cor- 
relator without the pion. On the bright side, the 
structure of the four-quark weak Hamiltonian re- 
quires a flavor change in the nucleon and therefore 
there are no disconnected diagrams to be com- 
puted on the lattice. 

The extraction of from N — > Nn re- 

quires an injection of energy at the PV weak ver- 
tex which can occur because the weak operator 
is inserted on one time-slice only. Therefore, we 
must include contributions from operators that 
are total derivatives, which usually vanish. Re- 
cently, chiral perturbation theory has been used 
to describe K — > ttt: with the kinematics ap- 
propriate for a lattice determination of the ma- 
trix elements of the relevant four-quark operators, 
m'jl** = m]^ u and m}^ tt — 2m 1 £ tt , including the 
necessary total derivative terms [^6|. In QCD, the 
LO Lagrange density describing PV interactions 
is given by 

/ 



r 



-h 



(i) 



N [ X 



X 3 n 



N (11) 



Vaa 



T, 



abd,fi 



(,(!) - 



h.c. 



while the Lagrange density at NLO is (keeping 
only nucleon operators) 



(NLO) 
wk 



= -2- iv ■ D N \ X 

A L 



xl 



N , (12) 



where i> M is the nucleon four-velocity. This is the 
leading contribution from a heavy baryon reduc- 
tion of iD^Nj^ [Xl - X%]N. Given baryon 



6 



number conservation, the total derivative gives a 
non-zero contribution from the energy and mo- 
mentum injected by the X£ — insertion. 
Working in the frame where the initial state nu- 
cleon (proton) is at rest, = (1, 0, 0, 0), the 
amplitude at NLO resulting from eq. (|12| ) and 
eq. (O) is 



•A-r, 



d 3 x C AI - 



\E)\p) 



-Vr, 



irNN 



E 
1 



(13) 



where E is the energy injected by the weak vertex. 
In order to produce an on-shcll wn: final state, the 
injected energy must exceed E > + M n — M p . 
Near threshold, where the final state neutron and 
pion are at rest and E — + M n — M p , the 
contribution from the total-derivative operator, 
and is formally dominant 



hp , scales as 



1/2 



over loop corrections and counterterms |lq ). 

3. Conclusions 

Lattice simulations will first make rigorous pre- 
dictions about observables by matching onto the 
appropriate EFT and using its explicit quark 
mass dependence to extrapolate from the lattice 
quark masses to those of nature. Therefore the 
EFT's must be known, which is the case in the 
meson sector and single baryon sector. However 
more work is required in the multi-nucleon sec- 
tor to be sure that the candidate theory J2?j (for 
a review see Ref. |2§|l ), in fact, is consistent and 
converges. PQ simulations represent the future of 
this field until fully unquenched calculations can 
be performed at the physical values of the light 
quark masses, and in this work we have made the 
small step of including baryons in PQQCD. We 
have constructed PQ^PT and shown that there 
are some interesting features beyond QCD. This 
field is just beginning and calculations beyond 
one-loop level are certainly required in order to 
understand the convergence properties of PQxPT 
and the uncertainties introduced in chiral extrap- 
olations. 

I am indebted to my colleagues, Silas Beane 
and Jiunn-Wei Chen, my collaborators on these 



works. 
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